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Abstract 

We consider quantum general relativity in three dimensions with a positive cos- 
mological constant. The Hartle-Hawking wave function is computed as a function of 
metric data at asymptotic future infinity. The analytic continuation from Euclidean 
Anti-de Sitter space provides a natural integration contour in the space of metrics, al- 
lowing us - with certain assumptions - to compute the wave function exactly, including 
both perturbative and non-perturbative effects. The resulting wave function is a non- 
normalizable function of the conformal structure of future infinity which is infinitely 
peaked at geometries where X + becomes infinitely inhomogeneous. We interpret this as 
a non-perturbative instability of de Sitter space in three dimensional Einstein gravity. 
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1 Introduction 

The application of quantum mechanics to cosmology is notoriously difficult. Even basic 
notions - such as unitary evolution on a Hilbert space - are difficult to reconcile with cos- 
mological features such as space-like singularities and inflationary expansion. Perhaps the 
most notable proposal in this regard is that of Hartle and Hawking p], who argued that 
the universe is described by a standard quantum mechanical wave function which can be 
computed via a Euclidean path integral. Without a complete understanding of quantum 
gravity however, this proposed wave function is difficult to compute precisely. The goal 
of the present paper is to investigate this wave function in 2+1 dimensions, where one can 
give a precise prescription for the computation of this path integral and understand its basic 
features. 

We will study 2+1 dimensional Einstein gravity with a positive cosmological constant: 



The equations of motion state that space-time is locally, though not necessarily globally, 
three dimensional de Sitter space. We wish to compute the wave function of the universe 
^(h), regarded as a functional of the topology and geometry (denoted collectively h) of a two 
dimensional spatial slice. Although the theory contains no local degrees of freedom the wave- 
function ty(h) still contains a great deal of interesting information about the global structure 
of space-time. In particular, it is possible to find classical solutions (such as quotients of 
dSs) with arbitrary topology and metric feona given spatial slice; *&(h) can be regarded as 
a wave function over this family of classical geometries. 

Schematically, the evolution of the wave function $?(h) should be determined by a 
Lorentzian path integral of the form 
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where the integral is over a suitable space of three dimensional Lorentzian manifolds M. The 
dependence on h enters through the boundary values of the three metric g. In addition, one 
must set some initial data which will fix the choice of state. For example, in perturbative 
quantum field theory around a fixed de Sitter background one obtains the usual Bunch-Davies 
wave function by demanding that the metric fluctuations are purely positive frequency across 
a past event horizon. At the non-perturbative level, however, it is unclear precisely how one 



should fix the initial data in (1.2). The proposal of Hartle-Hawking [JJ is that the boundary 
conditions are determined by a particular analytic continuation to Euclidean signature. More 
precisely, the wave function is given by a "no boundary" prescription: one integrates over 
Euclidean geometries whose only boundary is the spatial slice where the data h are fixed. 
According to the original proposal of pQ, each such geometry is weighted by the action of 
general relativity in Euclidean signature. 

The wave function *&(h) described above is also a central object of interest in the conjec- 
tured dS/CFT correspondence [2] (see also [31 H])Q As the spacelike slice is taken towards 
X + the volume of the spatial slice diverges and the phase of the wave function oscillates 
rapidly. In particular, near X + the wave function takes the form (see e.g. [31 E]) 

*(/0~e^W* ([h]) reQ . (1.3) 

Here S ct is a local functional of h which diverges near X + . The dependence of the wave 
function on non-local properties of h is captured by the complex function ([h]) . In 
particular, \l/ ren depends on the conformal class [h] of the spatial metric; all other information 



about the spatial metric is scaled away as one approaches X + . Equation (1.3) is the de Sitter 
analogue of the usual holographic renormalization procedure in Anti-de Sitter space (see 
[TO] for a review). In the dS/CFT correspondence \l/([/i]) rcn is identified with the partition 
function of a Euclidean conformal field theory. In the present case we have a purely metric 
theory of gravity, so consider only the dependence of \l/ ren on the conformal structure [h] of 
the geometry at X + . For theories of gravity containing additional bulk fields the partition 
function would depend on additional boundary data, and would be regarded as a generating 
function for correlation functions of the corresponding boundary operators. 

Our goal is to compute \£ ([/&]) directly in gravity, rather than via a conjectured boundary 
dual0 We will use an analytic continuation to Euclidean gravity. There are two distinct 
types of analytic continuation to Euclidean signature which could be used to do this. The 
first analytic continuation procedure is that originally used by Hartle and Hawking [1] (see 



also [llj). It relates the action (1.1) to that of general relativity in Euclidean signature with 



a positive cosmological constant. At the level of perturbation theory this can be understood 



1 See [5] for a more recent review and [HEIIH] for related approaches to holographic cosmology. 

2 In the following discussion ^([/i]) will refer to the renormalized wave function; we will omit the subscript 



in 11.3 
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t = 




(a) In the first type of analytic continuation, the (b) In the second type of analytic continuation, 
matching between de Sitter space and the sphere the matching between de Sitter space and Eu- 
happens on the equator at t = 0. clidean AdS (the hyperboloid) happens at X + . 

Figure 1 



as a complex coordinate transformation which takes the global de Sitter solution into the 
sphere. Explicitly, if t is the global de Sitter time coordinate, then under the Wick rotation 
t — > it: 



ds 2 



-dt 2 + cosh 2 tdfl 2 ) f dS (dt 2 + cos 2 tdQ 



:i-4) 



The surface of vanishing extrinsic curvature at t = of the Lorentzian geometry matches onto 
the t — slice of the Euclidean geometry. This is sketched in Figure (la). In perturbative 



quantum field theory, the resulting Euclidean path integral leads to the usual Bunch-Davies 
vacuum state used to compute fluctuations in inflationary cosmology. 

For our purposes, however, the standard Hartle-Hawking approach has two drawbacks. 
First, it is not known whether the full Euclidean gravity path integral with a positive cosmo- 
logical constant can be defined and computed precisely. Second, this procedure computes the 



properties of the state at finite time (t = in Figure (la)). In order to compute ^([/i]) one 
would need to evolve this wave function to the future boundary at t — > oo using a Lorentzian 



path integral of the form (1.2), which we likewise do not know how to compute precisely. It 



would be preferable to use an analytic continuation procedure which circumvents this inter- 
mediate step and allows us to compute \P ([/&]) directly via a Euclidean computation. Indeed, 
the holographic principle suggests that perhaps asymptotic quantities such as \& ([h]) can be 
well defined even if ^f(h) itself is not well defined at finite time. 

Both of these difficulties are resolved by the elegant analytic continuation procedure 
proposed by Maldacena [3]|^] The starting observation is that the analytic continuation 
z —7- —iz, where z is the planar time coordinate, along with ^Ads — ► ^ds takes de Sitter space 



i See also [TIH21I5] and references therein. 
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into Euclidean Anti-de Sitter space (i.e. the hyperboloid): 



d** = f MS r**'!** ) -+ 4 (^±^ ) . (i.5) 



Moreover, the usual Bunch-Davies boundary condition at the past horizon z — > — oo be- 
comes, upon analytic continuation, the condition that modes are normalizable in the interior 
of Euclidean AdS. Thus the wave function ^(h) is computed in terms of a path integral 
for Euclidean gravity with a negative cosmo logical constant. An important feature of this 
continuation is that the future boundary of the Lorentzian de Sitter geometry coincides with 
the asymptotic boundary of Euclidean AdS; they are both given by z = in (1.5). This is 
indicated in Figure (lb). This analytic continuation procedure computes the wave function 
at X + directly, circumventing the need for any Lorentzian computation. Under this ana- 
lytic continuation, equation (1.3) becomes precisely the usual relation between the partition 
function of bulk gravity in Anti-de Sitter space and that of a CFT, including the usual local 
counterterms. Thus this method is well adapted to the computation of the boundary wave 
function ^([/i]). 

The analytic continuations described above can be regarded as two attempts to define a 
path integral via a choice of contour of integration through the space of complex metrics. 
Without a more detailed understanding of this path integraQit is not clear which of these pre- 
scriptions is correct. At the level of perturbation theory around a fixed de Sitter background, 
however, both analytic continuation prescriptions appear to define the same state, at least 
for general relativity without additional degrees of freedom [3j. The question is whether 
they differ at the non-perturbative level. Indeed, it is not clear that the original Hartle- 
Hawking path integral can even be defined precisely beyond perturbation theory around a 
fixed classical background]^] Path integrals in Euclidean AdS are better understood. For 
three dimensional gravity with certain boundary conditions they can be computed exactly. 
Given that the Euclidean AdS results match those of the usual Hartle-Hawking contour at 
the perturbative level - and can be defined precisely at the non-perturbative level - we will 
therefore use the Euclidean AdS contour in the computation that follows. 

Our primary technical tool is the computation of the partition function of Euclidean AdS 
gravity by [TF] , which followed earlier work [18l [19], EOj [21] . The goal of the present paper is 
to reinterpret these results in the context of quantum cosmology]^] 

It is worth emphasizing that for certain theories of gravity the analytic continuation 
from de Sitter space to the hyperboloid is much more natural than that from de Sitter 



4 Such as that described in |T3l [33] for certain relatively simple theories. 
5 See [151 116] for an attempt to do so in three dimensional gravity. 

6 In [T7] the authors argued that the Euclidean AdS gravity path integral could not be interpreted as the 
computation of a canonical ensemble partition function in Lorentzian AdS. This conclusion was based on the 
traditional Hartle-Hawking analytic continuation prescription which relates Euclidean AdS to Lorentzian 
AdS. In the present paper we reinterpret the analysis of |17| as defining a state in Lorentzian de Sitter space, 
using Maldacena's analytic continuation prescription. Thus the conclusions are quite different. 
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space to the sphere. The isometry group of de Sitter space, 5*0(3, 1), is identical to that of 
Euclidean AdS. If gravity is formulated as a gauge theory - as, for example, in the Chern- 
Simons formulation of 3D gravity - SO (3, 1) appears as the group of local gauge symmetries 
[22] 123] . In a gauge theory, the usual Euclidean continuation prescription involves a complex 
change of coordinates but not a change of the gauge group. This indicates that the analytic 
continuation from dS to Euclidean AdS is in a sense the correct one, if one formulates gravity 
as a local gauge theory]^] The sphere, on the other hand, has isometry group SO (A) and is 
only obtained if one in addition changes the gauge group. 

In addition to three dimensional Einstein gravity, there is an second theory to which 
this observation applies: Vasiliev's higher spin theory of gravity in four dimensions [21]. In 
particular, Vasiliev theory with a positive cosmological constant has dS4 as its maximally 
symmetric Lorentzian solution, and Euclidean AdS4 as its maximally symmetric Euclidean 
solution. Much like the Chern-Simons description of three dimensional gravity, the theory is 
formulated purely in terms of connection variables. Indeed, recently there has been a great 
deal of progress on the formulation of the dS/CFT correspondence for higher spin theories 
of gravity [23 I2S1 [23 [2S], based in part on the analogy with the corresponding theories in 
AdS. The wave function of Vasiliev gravity with future de Sitter boundary was studied using 
the conjectured boundary gauge theory dual [29J, and shown to exhibit similar features to 
those we will discover below. 

This discussion highlights, however, an important drawback of the present approach, 
which is that it only applies to certain special theories of gravity. For generic matter coupled 
to gravity, the analytic continuation (1.5) makes the solution complex. Without a more 
detailed understanding of the path integral, it is not clear that such complex solutions 
should be regarded as genuine saddle point contributions to an appropriate path integral]^] 
The Maldacena contour should therefore be applied with caution. 

In the next sections we describe in more detail the exact computation under consideration, 
where the wave function is computed explicitly for geometries which asymptote to a torus 
at future infinity. We first describe the family of Lorentzian geometries on which this wave 
function is supported - essentially de Sitter analogues of the BTZ black hole - and then 
describe the explicit computation of the wave function. We will conclude that the wave 
function is non-normalizable, which is interpreted as a quantum mechanical instability of 
the de Sitter vacuum. 



7 We note, however, that the Chern-Simons formulation of 3D gravity does not appear to capture all of the 
desired features of quantum gravity, in particular black hole entropy |21j . Thus it is not clear how seriously 
one should take this argument. 

8 See [301 EU for a discussion of the contribution of such saddles. 
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2 The Classical Geometries 



The locally de Sitter solutions of Einstein gravity with a positive cosmological constant can 
have arbitrary topology at future infinity. In the maximally symmetric solution - global 
de Sitter space - X + has the topology of the sphere S 2 . For quotients of de Sitter space, 
however, the topology can be more interesting. The wave function ^([/i]) depends on the 
conformal structure of the metric h at future infinity; for the sphere this conformal structure 
is completely fixed. Thus on global de Sitter space there is not much to compute, at least 
in a theory of pure gravity without additional fields; the wave function on the sphere is 
just a number. To obtain interesting dynamical information about the wave function of 
the universe we should consider more complicated topologies for which the moduli space of 
conformal structures is non-trivial. 

In the semi-classical limit, the wave function \P ([/&]) should be regarded as a wave function 
over geometries which solve the equations of motion. We will first describe these classical 
geometries. We will consider in detail the case where future infinity has the topology of 
a torus T 2 , before commenting at the end of this section on the more general case. The 
conformal data of T 2 is encoded in a complex parameter r, so we denote the corresponding 
wave function \^(r, t)q 

The classical geometries for which future infinity has T 2 topology are quotients dS^j^L. 
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To see this simply, let us write the metric on de Sitter space in cylindrical coordinate^ 

d s 2 = -dt 2 + cosh 2 td9 2 + sinh 2 td(j) 2 . (2.1) 

In global de Sitter space ~ <fi + 2tt and 9 runs from — oo to oo. We can further identify 
9 ~ 9 + 27r to obtain a geometry which approaches a torus at future infinity t — > oo; this is 
an example of a quotient dSa/Z. The surface t = is a coordinate singularity in global de 
Sitter space, but in the quotient becomes a genuine cosmological singularity of Milne type. 
The wave function \I/(t, t) is supported on a family of geometries with big-bang singularities 
of Milne type. 

In fact there are an infinite number of such saddle point geometries dSs/Z which gen- 
eralize ( |2.1 ), which we will now construct]^] They will be labelled by elements of the coset 



Z\SX(2, Z). We start from the observation that Lorentzian dS 3 is the quotient SL(2, C) / SL(2, 
An element of this quotient can be decomposed as 

g = h(h*)- 1 , heSL(2,C). (2.2) 



9 We use this notation to emphasize that the wave function of Einstein gravity is neither a holomorphic 
function of r nor, necessarily, the square of one. In order to obtain a holomorphic function of t one would 
have to consider some version of topologically massive gravity or chiral gravity in de Sitter space, as discussed 
in e.g. [IS] . 

10 These coordinates are related by analytic continuation r — > ir, 9 —> iO to the usual coordinates on the 
static patch. 

11 These geometries have been considered before in the context of dS/CFT (see e.g. [32 ). 
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From here it is clear that g is invariant under h — >• ha with a G SX(2,R). The action of 
SL(2, C) on (7 and ft, 



Xh , Xg(X*y L , A G 5L(2, C) 



leaves the metric 



ds 2 = _j^ Tl{g -l dgg -l dg) 



(2.3) 
(2.4) 



invariant. 

The inflationary (i.e. planar) patch of de Sitter space is obtained by choosing coordinates 
for g 

w/z z — ww/z 
— l/z w/z 
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(2.5) 



with z G R and w G C. The metric is 



(2.6) 



In this language it is easy to construct the quotients dSa/Z. We make the identification 



g ~ 7fl ,( 7 



with 



7 



t = n + ir 2 



where (ti,T2) are real, 12 In the coordinates (2.5), the identification (2.7) implies 



(2.7) 
(2.8) 

(2.9) 



A more convenient set of coordinates is given by 

U\2 \ V2 



Z 



w 



t 2 - r 2 



exp(r 2 + ri6>) 



t 2 + T< 
t 2 - r 2 



2\ 1/2 



exp(— zr(0 + i#)) 



where now the action of (2.7) on the coordinates (t,<j>,9) is 

) + i^~u) + 27T~u) + 2"7TT . 



In these new coordinates the metric is 
ds 2 t 2 dt 2 



(t 2 -T?)(t 2 + T$) 



+ r 



TlT 2 
t 2 



de) + 



(t 2 + r 2 )(t 2 - r 2 ) 



t 2 



^de 2 



(2.10) 



(2.11) 



(2.12) 



12 The element 7 here is exactly the same as the one that defines Euclidean BTZ. See appendix [a] for 
further details on the relation between this construction and that of the BTZ black hole. 
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This describes a big bang/crunch geometry, with Milne singularity at t — T\. The conformal 
class of the spacelike boundary at X + (t — > oo) is a torus with conformal structure parameter 
r. We note that this geometry has no static patch. 



Equation (2.12) describes one possible geometry on which the wave function ^(r, r) is 



supported. In fact there are an infinite number of such classical geometries, each of which has 
a boundary torus with the same conformal structure r at future infinity. These geometries are 
related to (2.12) by modular transformation. To see this, note that the and 9 coordinates 



are not treated democratically in the geometry (2.12). In this construction we made a choice 
of contractible cycle; in the coordinates (2.12) the 9 circle shrinks to zero size at t = r 1; but 4> 



does not. The other geometries will correspond to other cycles which are made contractible 



at the Milne singularity. These geometries will be related to (2.12) by the action of the 



modular group SL(2, Z), which is the mapping class group of the torus. 

We can construct these geometries explicitly by performing a large diffeomorphism that 
acts non-trivially on the boundary. The modular transformation 



T — > 7T 



ar + b 
cr + d ' 
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d 



G SX(2,Z) 



is implemented by the diffeomorphism 



w — > 



w 



(2-13) 



(2.14) 



cr + d 

This provides a family of geometries labelled by SL(2, Z). We note that shifting r — > r + 1 



does not change the identifications (2.9). Thus distinct geometries should be labelled not 



by SL(2, Z) but by elements of the coset Z\SL(2, Z), where Z is the subgroup generated by 
T=(1\):t->t + 1. 

We emphasize that, although these new solutions are related by the coordinate transfor- 



mation (2.14), they are regarded as distinct contributions to the path integral. We adopt 



the definition that diffeomorphisms whose action on the boundary is non-trivial are physical, 
and that only those diffeomorphism that act trivially at future infinity are truly pure gauge. 
This definition was first discussed in the original proposal of the dS/CFT correspondence 
[2] . Further, following the construction of conserved charges in the dS / CFT correspondence 
(see e.g. [521 El]), it is straightforward to check that these modular transformation act 
non-trivially. 

An important observation is that these geometries are an incomplete classification of 
geometries whose boundary is a two-torus. Since we are not restricting the classical phase 
space to include only smooth Lorentzian solutions, there is a priori no reason to exclude more 
general quotients. For example, quotients of dS 3 by Z x Z m (with m an integer) satisfy our 
boundary conditions, and the additional action of Z m creates a conical singularity on top of 
the Milne singularity. However, as we will see below it is the quotients by Z which continue 
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naturally to the Euclidean BTZ black hole solutions via Maldacena's analytic continuation. 
Thus it is most natural to regard the wave function as supported on the quotients dS^fL. 
We will elaborate on this point in the next section. 

Given the above discussion, it is clear at least in principle how one can construct quotients 
of dS 3 for which X + has more complicated topology. In global de Sitter space the SL(2, C) 



isometries (2.3) act naturally on the boundary as conformal transformations. These are the 
usual Mobius transformations of the S 2 at future infinity. To obtain a non-trivial Riemann 
surface at I + one quotients de Sitter space by a discrete subgroup of SL(2, C), i.e. a Kleinian 
group. By uniformization, any compact Riemann surface with arbitrary conformal structure 



can be obtained by such a quotient of the Riemann sphere S . Equation (2.3) tells us how 
to extend this identification into the bulk to obtain a quotient of dS^. The conformal moduli 
of the Riemann surface then enter as the continuous parameters that label the choice of 
discrete subgroup. The wave function is then regarded as a function on these geometries. 

As in the torus case, there will be many locally de Sitter geometries whose boundary is 
given by the same Riemann surface. Indeed, one can act on each such geometry with the 
mapping class group (the higher genus analogue of the modular group jSX(2, Z)) to obtain 
an infinite family of geometries each of which asymptote to the same conformal geometry 
at future infinity. While it is difficult to write the metrics on these geometries explicitly, 
they will contain at the very least Milne singularities where cycles in the boundary geometry 
shrink to zero size. 

3 The Wave Function 

We now proceed to compute the wave function of the universe using Maldacena's ana- 
lytic continuation prescription. This will allow us to compute both perturbative and non- 
perturbative corrections to the usual Euclidean vacuum state. We will focus on the case 
where the boundary at future infinity is a two torus. As the conformal structure of the 
torus is modular (i.e. SX(2,Z)) invariant, \l/(r, f) will be likewise be an 5L(2,Z) invariant 
function of r. It can be viewed as a function on the fundamental domain of the action of 
SL(2, Z) on the upper half r-plane. 

According to Maldacena's analytic continuation prescription [3], a correlation function 
on X + of massless fields will equal precisely the answer obtained by analytic continuation 
from Euclidean AdS where we take the central charge 

c A ds = -> «c dS = I 3 - 1 ) 

to be pure imaginary. 

Under this analytic continuation Lorentzian gravity with a positive cosmological con- 
stant becomes Euclidean gravity with a negative cosmological constant. The wave function 
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\l/(r, f) becomes the partition function of three dimensional Euclidean AdS gravity with 
torus boundary conditions, which we denote Z(t,t). This partition function was computed 
in explicitly in [17] following earlier work [181 HH1 1201 IZE] • The result is 

(3.2) 



(c,d) = l 
c>0 



where 



and 



ar 



CT 



b 

~d ' 
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a b 
c d 



E SX(2,Z) 



Z ,i (r, r) 



q 



-fc+l/24M _ ? | 



|^(r)P 



(3.3) 



(3.4) 



Here 77 is the Dedekind eta function. The coupling k = CAds/24 is real for AdS gravity. 
Let us now describe the various terms in this expression and their interpretation in dS 



gravity. The summation in equation (3.2) is over all smooth Euclidean solutions which have 



a torus at their asymptotic boundary. These are the SX(2,Z) generalizations of the BTZ 
black hole, whose geometries are reviewed in appendix [Aj It is this sum which guarantees 
that the partition function Z(r, f) is a modular invariant function of r. 



The sum over (c,d) in (3.2) can be regarded a sum over the coset Z\5'L(2,Z), where Z 
is the subgroup generated byT = ( -J G SX(2,Z). Geometrically, T is a Dehn twist which 
leaves the bulk geometry unchanged. In the BTZ language, T is a modular transformation 
which does not change the mass and angular momentum of the black hole. In any case, the 
summand Zq i is invariant under T, so only by considering the sum over the coset Z\SX(2, Z) 
can one obtain a finite sum. We emphasize that this coset labels all smooth solutions to the 
equations of motion with a torus at the boundary. Thus all smooth saddle points which give 
non-perturbative contributions to the wave function are included]^] 

Z 0j i is the partition function of the thermal AdS saddle, including all perturbative (loop) 



corrections. The \q\ 2k term in (3.4) is the regularized classical action of Anti-de Sitter 



space, and the remaining terms describe a one loop correction to this action. The boundary 
conditions which define the de Sitter wave function in the far past are precisely those obtained 
by analytic continuation from the requirement of regularity in the interior of Euclidean 
AdS. Hence, around the vacuum solution the perturbative corrections to the wave function 
map precisely to the perturbative corrections to the classical partition function given by 



equation (3.4). In Euclidean AdS gravity these perturbative corrections can be computed 



in two different ways. The method of reference [17] involved a direct quantization of the 



13 In principle, it may be possible to alter the definition of our path integral to include non-smooth Euclidean 
saddle points. Without a precise definition of the path integral we have no way to know for sure. However, 
we do not expect that our qualitative result 
space - will be altered. 



the divergence of the wave function at boundaries of moduli 
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appropriate phase space of metrics. The corrections can also be computed using a more 



traditional one-loop determinant [35] • I n either case, the answer is (3.4). 

There is in principle a third method, which is to simply compute the corrections to the 
Hartle-Hawking state for perturbative gravitons directly in de Sitter space, as was done in 
[5]; this would presumably resemble closely the computation of [35J. We will focus on the 
Euclidean AdS construction, as only in this case are we able to include non-perturbative 
corrections. 

Indeed, the analytic continuation CAds _ > icds allows us to identify the wave function 
\&(t, f ) as the analytic continuation of the partition function Z(t, f) at the non-perturbative 
level. This is not merely a convenient computational tool. It provides a precise definition 



of the wave function via the Maldacena contour. In fact, the sum over geometries (3.2) can 



now be interpreted in Lorentzian signature as a sum over dSs/Z cosmologies given by (2.12) 



Each of the Lorentzian geometries described in the previous section will analytically continue 



to a smooth Einstein manifold in Euclidean AdS. Thus each of the terms in the sum (3.2) 
can be interpreted directly as a quantum mechanical amplitude associated to a particular 
cosmology. 

In principle one should worry that the form of these non-perturbative corrections is subtle, 
in that the answer differs depending on whether we analytically continue in CAds before or 
after performing the sum over geometries. From the explicit results below, however, this 
does not seem to be the case. We will therefore perform the sum first and then analytically 
continue, using the results of [T7j. 



We first consider the analytic continuation of the leading perturbative piece (3.4) when 
c Ads - > icds- This means that k = icds/^ is purely imaginary. In the usual perturbative 
regime where r = t\ + ir 2 — > ioo the partition function Z i diverges like 



J 0,l 



g|- 2fe (lH ) = e 4 ^ 2 (l + ---) , (3.5) 



when k is real. This becomes a pure phase when k is imaginary. Thus the classical contri- 
bution of Z 0) i to the wave function ^(r, r) is a pure phase at r — > ioo. However, the next 
to leading piece also diverges at r 2 -> 0, because at this point 

M _ p -2ttt 2 \2 

Thus even though the perturbative piece is finite at the "usual" pole at r = ioo, it has a 
divergence at the modular image of this pole r = 0. In Euclidean AdS (with real central 
charge) this divergence is of no concern; it is much smaller than a stronger divergence which 
will occur due to one of the other saddles (the Euclidean BTZ saddle) at r = 0. In the 
Euclidean AdS picture this divergence can be interpreted as a one loop renormalization of 
the vacuum energy of the ground state. In de Sitter space, however, this divergence is crucial. 
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Unlike the classical divergence (3.5), there are no factors of k in (3.6). Thus while the usual 
tree level divergence (3.5) becomes a phase under analytic continuation the one loop piece 
does not. 

In fact, because the Dedekind eta function is (almost) modular invariant we can argue 
that Z 0; i will have a pole at every rational point on the real r axis. This is because in (3.4) 
the |(z| 1,/12 factor cancels the pole in the eta function at r = zoo, but not at any of the modular 
images of zoo. Based on this, it is clear what should happen when we sum over geometries; we 
will find that the wave function has a pole at the cusp r = zoo, because of analogous one-loop 
terms which appear in the expansion around other (modular transformed) geometries. 

In order to confirm this expectation, we must discuss an additional subtlety. The partition 
function Z(r,f) is formally divergent due to the sum over Z\5X(2, Z). However, this sum 
can be regulated in a natural modular invariant way. The result is complicated but can be 
evaluated explicitly in some cases. To write the answer in a simple way, let us note that 
Z(t,t) is periodic under r — > r + 1. So we can Fourier expand in t\. Let us just consider 
the T\ independent piece. Then completely explicit expressions can be found for the wave 
function following the derivations in [T7] ; see appendix B for explicit formulas. In particular 
if we expand near r = zoo, the terms independent of T\ are given by 



* = £o,i 



1 oo 

- 2 J2c n mr 2 - n , (3.7) 



n=0 



where Zq^ is the analytic continuation of (3.4) to imaginary k, and c n are constants which 
can be computed explicitly. They are polynomials in k of order n. For example 

(7r 3 -67r)(ll + 24fc) 
Co = -6, Cl = m . (3.8) 

As anticipated above, the wave function has a divergence at r — > zoo of the form 

^(r, f) ~ — 6 exp | — Imr| as Imr — > oo . (3-9) 

This divergence is precisely the one anticipated above; the factor of —6 is due to the regular- 
ization of the sum over modular images. We emphasize that Zq^, which includes the classical 
action of the usual (thermal AdS) saddle along with its perturbative corrections, is finite. 
The divergence comes from non-trivial saddles, and in fact from the one-loop contribution 
to the action of these saddles. 

We conclude that the wave-function of the universe has a peak at very "asymmetric" 



universes where Imr is large, This peak renders the wave function non-normalizable 
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14 Recall that Imr is essentially a measure of the inhomogeneity of the torus. A rectangular torus with 
cycles of length 1 and has t = 1T2 . 

la This agrees with the computation of |15j . in which the wave function of C1S3 was considered using a more 
standard Hartle-Hawking path integral over Euclidean geometries with positive curvature. The advantage 
of the present computation is that we now have a physical interpretation of this divergence. 
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Thus the probability density is infinitely peaked at infinitely inhomogeneous universes. We 
note that this peak is exponential in Imr as Imr — > oo, thus the wave-function is non- 
normalizable even with standard SL(2,M.) invariant measure on the upper half r plane. 
It is also interesting that the strength of the peak (the coefficient in the exponential) is 
independent of the de Sitter radius. This is because it is essentially a quantum mechanical 
effect - due to a one loop effect, albeit around a "non-perturbative" saddle - rather than a 
classical one. 

Our interpretation is that three dimensional de Sitter gravity is quantum mechanically 
unstable. The wave function resembles that of a quantum mechanical system with a potential 
which is unbounded below. It is tempting to speculate that this mechanism is a simple 
version of the more complicated instabilities of de Sitter vacua observed in string theory 

In order to make this analogy precise, however, one would 
want to determine the endpoint of the instability. In the present case the answer to this is 
unclear. It might be necessary to add local degrees of freedom in order for the theory to have 
a sensible true vacuum ground state. Vasiliev's theory of higher spin gravity is a concrete 
setup to explore these issues [25] . 

Of course, another natural direction is to explore to what extent these results can be gen- 
eralized within dS3 gravity. Consider, for example, the wave function with more complicated 
asymptotics. When the boundary has higher genus, the corresponding Euclidean geometries 
are quotients of Euclidean AdS by Kleinian groups, and the resulting partition functions are 
more difficult to write down (see e.g. [101 SI])- It seems likely that in this case that the wave 
function will be infinitely peaked at the boundaries of moduli space where the geometry of 
X + is infinitely inhomogeneous. It would be interesting to investigate this further. 
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A Euclidean BTZ 

In this appendix we review the construction of the Euclidean BTZ black hole as the quotient 
H3/Z. This should be compared with the analogous construction in section [2] of dS^/T,. 

Our instability might even be related to the instabilities of de Sitter discussed in [37 ) l38l 139] . 



constructions (see e.g. [36J) 
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Euclidean AdS3 (i.e. H 3 ) is the quotient SL(2,C)/SU(2). An element of this quotient can 
be decomposed as 

g = hh ] , heSL(2,C)- (A.l) 

From here it is clear that g is invariant under h — >■ ha with a G SU(2). SL(2, C) acts on g 
and h as 

h^Xh , g^ XgX 1 , A G SL(2, C) . (A.2) 

The metric is 

ds^^-Tr^- 1 ^^ 1 ^) . (A.3) 

Note the sign difference with respect to de Sitter; they are related by the analytic continuation 
^Ads - > i^ds- In Poincare coordinates the element g is 



p + ww/p w/p 
w/p l/p 



and the metric is 



— Trig 1 dg g 1 dg) = — (dp 2 + dwdw) 
2 p l 



The identification that defines the BTZ solution is 



with 



7 



-IITT Q 



, T = Ti + IT 2 • 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



e lw '- 

This is the same element used for de Sitter, i.e. TAds = T ds- In the coordinates (|2.5l), the 



action of the identification (2.7) implies 



p ~ e - iw{T - f) p , w ~ e~ 2TiT w 



(A.l 



A more convenient set of coordinates is given by 

12 \ 1/2 

/ 1 ' 

P 



r z _|_ T z 



2 , _ 2 1 exp(r 2 + r x t) , 



r 2 + T 2 



2\ 1/2 



exp(-ir(0 + it)) , 



where 

w) = + it~w) + 27r~W + 27TT , 

in terms of which the metric reads 



r 2 dr 2 



ds 2 



^ids (r^-rD^ + r!) 



+ r 



dt) + 



(r 2 - r 2 )(r 2 + t 2 ) 



dt 2 



(A.9) 
(A.10) 
(A.ll) 



The Euclidean BTZ black hole geometry described above is related to the dS3 cosmology 



(2.12) by the analytic continuation ^Ads — ► ^ds an d r — >■ ir. 
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B Euclidean AdS3 path integral and its analytic con- 
tinuation 



In this appendix we collect the results for the Euclidean AdS3 path integral computed in 



|17j . The sum is given by (3.2) which is divergent due to an infinite set of saddles with 
arbitrarily small volume in Planck units. This divergence can be regulated in an SL(2, Z) 
invariant manner, giving 



Z(r,f) 



1 



\v(r)\ 2 



E(2k - ^,0) + E(2k + 2-^,0) 
-E(2k + 1-^,1) -£(2* + 1-^,-1) 



where CAds = 24fc. We have defined 

E(k, h) = e 2 ^ T2+l » Tl) + e- 2winT1 E n {K, fi, r 2 ) , 

n 

with r = T\ + ir 2 . The coefficients E n for n^O are given by 

oo 

where 



m=0 



and 



c=l 



de(z/cZ)* 



XmM = E c~ 2{m+1/2) S(-n, ^ c) , S(a, b; c) = £ exp ( ^(ad + bd' 1 ] 



For n = in (B.2) we have 



m 
2 ' 



m=0 



where 



w m (K, 0) 



2 m 7T m+1/2 C(2 



m 



mr(m + l/2)C(2m + 1) 
2 vr m + 1 /2 

Wm(K,±1) = mr(m + l/2)C(2m + l) Tm(K) 
and T m (x) is the Chebyshev polynomial of the first kind. 



(B.l) 



(B.2) 



W«,^r a ) = / ^e^(l + Z 2 )- m " 1/2 («- W) m (B.4) 



(B.5) 



(B.6) 



(B.7) 
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Figure 2: The norm |\I/(t)| of the wave function (B.8) for purely imaginary values of r (n = 0). 
In order to obtain stable numerical results it was necessary to include the ranges n = — 3 . . . 3 and 
m = ... 30 in equations (B.6), (B.3). 



The coefficients denning the regulated sum (B.l) are well defined for any complex k, so 



the regularization technique used in [T7] is sufficient to define the partition function in de 
Sitter space by analytic continuation in k. Taking CAds — > ic<is, the wave function is 



1 



\v(t)[< 



E(2ik-^,0)+E(2ik + 2-^,0) 
-E(2ik + 1 - ^ 1) - E{2ik + 1 - ^, -1) 



(B.8) 



To illustrate the shape of the wave function, in figure [2] we plot the sum of first few hundred 
terms in the series expansion of (B.8) for Cds = 2 (k = 1/12). Even though the coupling is 



small, the divergent behavior of ^(r) is already evident from explicit evaluation of the first 
few terms in the sum. 
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